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Abstract. 

In this article we review recent generalisations of the central limit theorem for the sum of 

specially correlated (or q-independenf) variables, focusing on q > 1. Specifically, this kind of 

l 

correlation turns the probability density function c £ q {X) = srf q 1 + (q — 1 ) f5 q (X — jX q ) 2 q , which 

emerges upon maximisation of the entropy S q = k (1 — / [p (X)] q dX) / (1 — q), into an attractor in 
probability space. Moreover, we also discuss a ^-generalisation of a-stable Levy distributions which 
can as well be stable for this special kind of correlation. Within this context, we verify the emergence 
of a triplet of entropic indices which relate the form of the attractor, the correlation, and the scaling 
rate, similar to the ^-triplet that connects the entropic indices characterising the sensitivity to initial 
conditions, the stationary state, and relaxation to the stationary state in anomalous systems. 
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INTRODUCTION 

In our previous article [1], we have verified the standard central limit theorem and 
its Levy-Gnedenko extension for the case of independent and identically distributed 
random variables associated with a ^-Gaussian distribution, 

%{x) = ^ q e q ^ h) \ (1) 

i 

with e* = [1 + (1 — q) x] (if 1 + (1 — q)x > 0, and zero otherwise) [e\ = e x ) . Dis- 
tribution (1) optimises the continuous version of the nonadditive entropy [2] S q = 

( w a\ 

k[l— Ei ? f)/( 1— where q G 91. In this article we review recent generalisations of 

the central limit theorem which have been formulated within nonextensive statistical me- 
chanical concepts. After numerical indications suggesting the existence of a generalisa- 
tion, within nonextensive statistical mechanics, of the central limit theorem [3,4] 1 , such 



The models introduced in references [3] and [4] have recently been analytically solved in Ref. [5]. It 
was verified that, for these two cases, the limiting distributions are not in fact g-Gaussians, but are instead 
other distributions. These distributions are, however, so close to g-Gaussians "that it becomes extremely 



generalisation has indeed been proved for variables with finite or infinite ^-generalized 
second-order moment [6, 7]. 



CENTRAL LIMIT THEOREMS FOR ^-INDEPENDENT 

VARIABLES 



The g-Gaussian case 

As we have reviewed and illustrated, in the case of convolution of independent random 
variables (including random variables associated with %(X) distributions), there are 
only two stable forms in probability space, namely, the Gaussian and the a-stable Levy 
probability density functions. We have also referred in Part I [1] to other versions of 
the central limit theorem that are available in the literature. In what follows we discuss 
the addition of random variables that are correlated in such a special way that a new 
algebra, the ^-algebra [8, 9], is necessary for a suitable analysis. In this context, it has 
been introduced in Ref . [6] a non-linear integral transform, the g-Fourier Transform, 
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,f(X)dX. 
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Applying the definition of ^-product, x® q y = [x l ~ q +y l ~ q — l] in Eq. (2), it 
is possible to write & q [/] (k) without the explicit use of the g-product, & q [/] (k) = 



r./POexp, 



ikX 



{/(*)} 



l-q 



dX. As an application of & q [f] (k), it is provable that, 



& q [%(X)] (*) = W 



2(9-1) 



3-q 



(3) 



where C q = \f$ j srf q , and j8 = SB with si and & as defined in Eq. (9) of Part I [1]. 

Another distribution for which it is simple to obtain its g-Fourier Transform, is 
the uniform distribution, ty/ (X), % (X) = ^ (—a < X < a, a > 0) 2 . Its ^-Fourier 

ak 



Transform is, & q [<% (X)] (k) = 

e c7-generalisation of sir 
Introducing a function, 



ak 



sm q 



the c7-generalisation of sin (x) [10], and q = 2 — i. 



q{2a) 



(«-l)/4 



where sin ? (jc) represents 



v(u) = 



l + u 



(u<3), 



(4) 



difficult to distinguish the true curve from its q-Gaussian approximant" [5]. These interesting results show 
that the probabilistic correlations included in these two specific models do not correspond exactly to q- 
independence, but are only very close to it instead. 

2 In this case, as well as for all distributions with compact support, integration must be done over that 
support. Otherwise the integral does not converge. 



whose inverse is, 



3m -1 
l + u 



(u>-v, 



(5) 



and assuming q\=v (q) and = v 1 (q), it is possible to rewrite Eq. (3) as follows: 



& q [%(X)] (k)=cxp qi [-P'k 2 ] 



(6) 



and & q _, (X)} (k) = exp q k 

A 

Equation (7) can be rewritten as 

i 



, where 
3 — q 



" %p-ic 2 q [q - iy 



(7) 
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2(9-1) 



q< 2. 



(8) 



We might consider the case q — 1 in Eq. (8). In this situation, the Fourier Transform 
of <S\ (X) has the same functional form, <S\ (k) = exp [— . For Gaussian functions 
like these two, the width (and the inflexion point in linear-linear scale) are related to 
/3 (actually ~^)- Hence, Eq. (8) reflects a relation between uncertainties in real and 

reciprocal spaces. In the general case, relation (8) is a sort of ^-analogue of the quantum 
mechanical uncertainty principle by Werner Heisenberg [1 1]. In Fig. 1 we represent 
K(q) for values of q between —5 and 2. 




FIGURE 1. Representation of Eq. (8) for values of q between —5 and 2. For q = 1 K(q) = 1/4; 
lim^-ooA^g) = 0. We are presently focusing on q > 1. 



Equation (6) has permitted to verify the mapping, through & q , 



■J. 



9 ' 



qi=v(q), l<q<3, 
q-l=v~ l (q), l<q, 



(9) 



9— 1 



and to prove the existence of the inverse q-Fourier Transform, & q 1 , c S qx — > 

(qi—v(q) : l<q<3), and ^ — > i (q^i =v^ 1 (q) , 1 < q). It is worth to men- 

1 

tion that and q_\ fulfil the dual relation q_\ + — = 2, that has also appeared in the 
context of phase space self-invariant occupancy [12], i.e., such as that all marginal prob- 
abilities of the system composed by N equal and distinguishable subsystems coincide or 
asymptotically approach (for N — > °°) the joint probabilities of the (N — l)-system. If we 
consider a sequence of applications of v (q), q n = v n (q), it can be seen that 

2 2 

■+n, n = 0,±1,±2, . . . and q = vq (qo) . (10) 



1 -q n 1 

If q = 1, then q n = 1 for all n. Otherwise, i.e. if q ^ 1, in the limit n — > ±oo, q n — > 1. This 
result can be interpreted in the following way: the simple successive application of the 
q-Fourier Transform, 

^ = m=l,2,...and n = 0,±1,±2, . . . (11) 

on a distribution % leads to the Gaussian form, lim &™ VS q \ = <S . 

Let us now present the scheme within which the functional form Eq. (1), i.e., 5f 9 (X), 
turns out to be stable. 

Two random variables, X and Z, are said q-independent if 

& q „ [& (X + Z)} (k) = ^ [p (X)] (k) ® 9 [p f (Z)] (k) , (12) 

where q* = q-\, (X + Z), p (X), and p' (Z) are the probability density functions for 
X + Z, X and Z, respectively. As an exemple of ^-independency we refer two variables, 
X and Z, both associated with a g*-Gaussian probability density function, <g q% (X) and 
^ (Z), with j8x and /3z, respectively. If the variables are ^-independent then the sum 



i_ 



X + Z is also associated with a q* -Gaussian whose width is 8 = f 3 q | 

Where ^ = 

Consider now the variable, Y = ~^-~^j^7^ iV/Jg * , where X,-, z e ( 1 , N), are identically 
distributed random variables. 

Let Xi, . . . ,Xjv be a sequence of ^-independent identically distributed random vari- 
ables with finite g*-mean, fi qt , and finite second (2g* — 1) -order moment, of^j. Under 
these conditions, with 

Av (<?*) = (V^^-i^-i) 2 ^ , (13) 

^(y) is said to be ^-convergent to a g_i -normal distribution as Af — > °°, with 3f u = 
I-oc [f (X)] u dX. Analytically, this can be written as, 

l^(Y)] (k) = %(k) = & q _ x (X)} (k) . (14) 



The proof of Eq. (14) [6] follows along the lines of the standard central limit theorem 
where, instead of & [/] (k) and the usual product, we use [f] (k) with 

[/] (*) = ! + W q * &q.* ~ f ^-i fc 2 + ^ (^ 2 ) , (* - 0) , (15) 

and the g-product. Thus, the convolution of ^-independent random variables Y is written 

as, 

^ [& (Y)} (k) = [p (X)] (k) ® q ...® q & qt [p (X)] (k) , (N factors), (16) 

where 1 < q < 2 and consequently 1 < q* < | . 

Using Eq. (16) in Eq. (15), together with properties of g-logarithm when k — > , 

In, & q _ t \& (Y)] (k)=Nln q (l-3?$ + (g) ) , i.e., 



N 



limJ^i [^(7)](*) = exp 9 



(17) 



In other words, Y is q-i-convergent to the random variable Z whose q-\-Fourier 
Transform is given by Eq. (17). Hence, according to the mapping relations, the explicit 
form of the corresponding q-\ -Gaussian, <3 q i (X), yields 

% ~' m = ^ (2VTT^f^ r "-' hft -' X2] ' <18) 



T- 



where j3 9 _, = (3 — q-\) / {Aq^ q ^ 1 ^ j . Obviously, when we ^-convolute q*- 

Gaussians, the resulting probability density function is a g*-Gaussian, i.e., ^ (X) is a 
stable attractor for q-independent random variables upon the condition stated above, 
in the same way Sf (X) and L a (X) are the stable attractors for the sum of independent 
random variables. 

Concerning the exponent of Eq. (13), it is easy to verify that (2 — g*) -1 = 

v(v(g_i)) = q\. Defining 6 = (2 — ^ l ) 1 , as the scaling rate exponent, we have 
8 = q\. We are then allowed to define, for generic n (following Eq. (10)), a ^-triplet, 
{qatt,qcomq S cal}> which relates entropic indices for the attractor, q att = q n -\, the 
correlation, q corr = q n , and the scaling, q scal = q n+ \. 

Consider now the variable Y' = YDpj(q*), that is composed by the sum of /V q- 
independent random variables all following a g*-Gaussian with the same j8^. Using 
the associative property of the g-product in Eq. (16), we are able to obtain the scaling 
relation, 

%.(?)= NP^, (19) 

in Fourier space (see Eq. (7)). For the standard (q* = 1) and generalised (2g* — 1)- 
variances, we have obtained the relations 

a 2 (Y)=N^a 2 , and a^_ x (Y) = N^o^, (20) 



with a 2 = [JB 9 , (5 -3q*)] 1 and a^_ x = [jB 9 . (1 + #*)] ^ When # = 1, Eq. (20) turns 
out into the well known relation for the sum of independent variables, o 2 (Y) = No 2 . 
The way on which variance scales upon addition is a standard tool to evaluate the 
character of a time series whose elements X ; have a variance o 2 . It is well known that 
the nature of a signal is characterised by its Hurst exponent, H, = N 2H o 2 where 
Ojy represents the variance of a new variable obtained from the sum of N elements 
of the time series with variance o 2 . The series is considered as anti-persistent if < 
H < 1/2, Brownian if H = 1/2, and persistent if 1/2 < H < 1 [13]. By comparing the 
Hurst exponent definition with Eq. (20) we verify that a connection can be established. 
Expressly, and from Eq. (20), we conjecture that ^-independent stochastic signals should 
respect the following relation H = [2 (2 — q*)]^ 1 . 



Verification of the q- generalised CLT 

We verify here that correlations of the form of Eq. (12) yield a ^-Gaussian attractor. 

Let us start with the case of g*-Gaussians with q* = 3/2 and /3 = 1 which are q- 
independent. By direct evaluation of its q* -Fourier Transform, J? qt [@ qt (X)] (k), and 
taking into account the Cauchy principal value, we obtain 

3 

(21) 

which corresponds exactly to a ^5 (k) function with the same q and /3 as indicated by 

Eq. (7). Using Eq. (21) in Eq. (12), and the mapping relations (4) and (5), we have ob- 
tained the convolution of two ^3 (X) distributions which is also a |-Gaussian. From the 

latter, and applying the associative property of the ^-product, we have calculated the 
convolution of N = 2,4, 8, 16 [q = |) -Gaussian distributions. Contrarily to what hap- 
pens for (q = 1) -independent variables, the resulting distribution is always a [q = |)- 
Gaussian that will never converge to a Gaussian, even when iV — > °°, and despite the 

finiteness of the ^3 (X) standard deviation. On the panels of Fig. 2 we depict the be- 

1 

haviour we have just described. 

Another illustration is presented in Fig. 3 for the case of the sum of (q — ^)- 
independent random variables associated with a ^9 (X) distribution with /3 = 1. As we 

have verified when the random variables have the same probability density function 
but are (q = 1) -independent instead, the convolution leads to a a-stable distribution, 
L a (y^oo), with a = |. In the case of [q = |) -independence, the limiting (stable) dis- 
tribution is the ^-Gaussian in accordance with the ^-generalised central limit theorem. 
Consistently, (N) = N 5 and j3; (N) = N fi' qt (1). 
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FIGURE 2. Upper panels: Probability distribution P(Y N ) vs. Y N , with Y N = Y^ =l Xi , X t being (#=§)- 
independent random variables associated with a £f 3 (X) distribution with j3 = 1 (fe/f), and the respective 

(<7 = |)-Fourier Transform, P (k), vs. k (right). Middle panels: Same as above, in ln3 -squared scale (left), 

and In 5 -squared scale (right). The straight lines indicate that P (Yn) and P (k) are g-Gaussians with q = | 

and <7 = | respectively. Their slopes are J3~^ 3 ; 2 (-AT) for left panel curves and (AT) for right panel 

curves. Lower panels: P~ =3 / 2 (^0 V,S- wr, i cn i s a straight line with slope 1 (left); =3 y 2 (N) vs. N 



which is also a straight line but with slope 



3-« 



0.088844... (ng/zf). 



?*=3/2 



The (<?, a) -stable distribution 

Within nonextensive statistical mechanics, the a-stable Levy distribution has also 

been generalised [7]. A distribution f(x) has been considered whose asymptotic form, 

i+g 

|X| — > 00, corresponds to /(X) ~ C \X\ l + a (i- 1 ) (for g = 1, f(x) recovers the already 
mentioned a-stable Levy distribution). Within this new class, we can define, e.g., a 
g-Cauchy distribution when a = 1, which leads to the classical Cauchy distribution 
for q = I. Just as the usual Levy distribution, this q-generalisation of the a-stable 
distribution, ^ q ,a (X), is defined by its q-Fourier Transform 

& q [Sf qta (X)] (k) = C'exp 9l \k\ a ] , (22) 
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FIGURE 3. Upper panels: Probability distributions P (Y N ) vs. Y N , with Y N = £f = t Xj , X t being (q = | ) - 
independent random variables associated with a (X) distribution with j3 = 1 (left), and the respective 

= |)-Fourier Transform, P (k), vs. k (right). Middle panels: Same as above, in Ins -squared scale (left), 

NUUU1UU NLiilU I / }. 1 JIL" MU1IUIU 11H.I11.IIC I1U1I r \ I Y 1 dllU r \K i ill L" f /~ v_Ji}U NNlill 1*. Willi (/ 



and ln7 -squared scale (right). The straight lines indicate that P (Y/v) and P (k) are g-Gaussians with q — j 
and q=\ respectively. Their slopes are P q=9 i 5 (N) for left panel curves and =9 ^ 5 (N) for right panel 
curves. Lower panels: —9/5 (^0 vs. N 5 , which is a straight line with slope 1 (left); =9/ , 5 (N) vs. N, 



which is also a straight line, but with slope 
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0.030995 . . . (right). 



where, as stated previously, q\ = =v(q). Distribution J£q,a P0 presents an infinite 

(2q — l)-variance, when l<g<2, 0<oc<2, and a < j^—. From this point on we 

denote, ^ q ,a (z) = Aexp q [— c \z\ a ] ■ Parameter A is related to the normalisation of the 
distribution and c basically controls its width. 

Along the lines of the previous generalisation, it has been shown that the sum of q+- 
independent random variables, all having the same distribution f(x), leads to Eq. (22). 
Since f(x) is stable, i.e., after an appropriate scaling, the sum of X variables has the 
same form for the probability distribution. In other words, f(x) is an attractor in the 
probability space, being a Jz? q , a (X). In addition, distribution f(X) is asymptotically 
equivalent to ( S qL (X) = ^ qL) 2 (X), where qL = (2qa — a + 3) / (a + 1). In this way, we 
can say that the following mapping has been asymptotically introduced, 2 — > ^ q ,a ■ 



In the second paper by the same authors [7], an extension for all 1 < q < 2, < a < 2 
has been introduced. This extension is based on the asymptotic equivalence between 

f(X) and the (q, a) -distribution, %, a {X) ~ \X\~^ (\X\ -> oo) 3 , and the fact that, 
for < a < 2, and arbitrary q\, there exists an index q2, together with a one-to-one 
mapping, such that & a (X) — > ^fq 2 ,2 (X). Using this mapping, it has been obtained 

the extension of Eq. (9) for ^ a (X), which is written in the following way 

%,a^%'.a, Kq<2,0<a<2, (23) 

where 

a , a , * a-2(l-o) 
_ =1 + _ and? . = (24) 

(a) standing for asymptotic behaviour. Moreover, it has been shown that this inverse 
g-Fourier transform exists, i.e., 

%>,a ^ %,a- (25) 

The consecutive application of n ^-Fourier transforms, ^ q * a , each one following 
Eq. (24), leads to, 

a a + n, n = 0,±l,±2,... (26) 



1 - q n l-q 

Hence, following the previous scheme, it has also been proved [7] that, if we consider the 
sum Y = Xi+) ^^ +Xn °f symmetric variables X t mutually q\ -independent, and having 

_ a+l 1 

a probability density function f(X) ~ |X| (|X| — >• oo, Dtv (g) °< N a( - 2 -^ ), then, 

when TV — > °°, ^ (J 7 )] (fc) = exp^ [— |fc| a ] , a (gi, oc) -stable distribution. Bearing in 
mind Eq. (24) and Eq. (25), we have that Y is convergent to a ^(^a) PO distribution with 

2(l-g)-«(l + g) 
* 2(l-*)-a(3-*)' U/) 

Regarding scaling, D^{q), assumed for variable Y = Xi+X ^^~ Xn ^ we can define the 

scaling rate exponent for 5f 9;Ct (X) (asymptotically equal to f(X)) as, 5 = [a (2 — g)] -1 . 

Summarising, we have two different approaches to the attractor of the sum of ran- 
dom variables following a (q, a) -stable distribution. These two approaches can be 
understood as the existence of a crossover between two regimes in the attractor for 
(q, Ot)-stable distributions. The first regime corresponds to the intermediate region of 
the attractor in which it is asymptotically equal to a ^ ;Ce (X), whose tail exponent is 



It should be emphasised that Sf 9)Ct (X) is not a (q, a)-stable distribution since it does not respect Eq. (22). 



[2 (1 — q) — a (3 — q)\ / [2 (1 — q)\. When a = 2, this exponent coincides with the expo- 
nent of the attractor for the convolution of g-Gaussians, 2/(q—l). The second regime, 
which tends to the a-stable Levy distribution when q=l, represents the behaviour for 
large |X| where the attractor, S£ qM (X), is asymptotically equivalent to a &q L ,2 PO distri- 
bution which has a tail exponent equal to (oc+ 1) / (1 + (Xq — a). When q = 1 this slope 
coincides with the exponent for a a-stable Levy distribution, a + 1. Yet in this regime, 
and q = 1, we verify qi = which coincides with the relation obtained in reference 
[14]. For a — > 0, we have qt — > 3, upper limit for normalisation of <3 q (X), and when 
a — > 2, <£l — > |, z'.e., the upper limit for finiteness of the second-order moment, of, of 
%{X). 

The former analysis is evocative of Fig. 5 of Ref. [1]. In other words, as a — > 2, 
=^,a (X) approaches the ^-Gaussian, ^(X) and, at some critical value X c , it changes 
its behaviour to a power-law decay with exponent (a+l) / (1 + aq — a). With this 
picture in mind, we sketch in Fig. 4 the attractor for the case q\ = 2 = v(g = 4) and 
values of a approaching a = 2. We can verify that, for all 1 < q < 2, the inequalities 
_2_ y 2(i-q)-a(3-q) > _g^L_ hold. Our sketch might be corroborated in the near 

future, as soon as the form of the inverse g-Fourier Transform becomes analytically 
available. A summary of the whole situation is presented in Table 1. 
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FIGURE 4. Outline of (q, a)-stable distributions for the case in which the correlation is given by =2. 
As a of Levy distributions approaches 2, the distributions (q, a)-stable becomes more and more similar 
to a G| (X) with an exponent [2(1 — q) — a (3 — q)} / [2(1 — q)]. However, since a^2, for some critical 

value X* , the distribution changes to another regime with a tail exponent (a + l)/(l + a^ — a). 



FINAL REMARKS 

In this article we have reviewed the generalised central limit theorems presented in 
[6, 7]. These theorems are based on nonextensive statistical mechanics and they address 
the sum of q\ -correlated random variables (with q\ = whose attractor scales as 
jy-\/\a(2-q)]_ introducing an appropriate nonlinear generalisation for the Fourier Trans- 
form, it has been possible to verify the existence of a new attractive subspace in prob- 
ability space, namely £f 9;Ce (X), which contains the Gaussian, <S (X) = ^2 (X), and (in 
an asymptotic way) a-stable Levy distributions, Jz?i ; a (X) = [^i,a (&)] (X). These 
special attractors are both related to the sum of independent variables. Within £f 9;Ce (X), 



TABLE 1. Resume of the main results presented in the article: Central limit theorems which present a N 

l/[a(2-?)]. scaled attractor F(X) 

for the sum of N — > °° ^-independent identical random variables with symmetric distribution f(X); q\ = P^-; we focus on q > 1. The term 
intermediate must be read as an infinity, however not so large as the infinity associated with the distant regime. For q ^ 1 and a = 2, when 
the random variables are associated with a ^-Gaussian, we verify the scaling relation, j5'(N) = Nfi', where j3' is the (inverse) width for the 
g-Fourier Transform. 





q = 1 [independent] 


q 1 [globally correlated] 


02q-l < °° 

(a = 2) 


F(X)=&(X) 
with same a\ of / (X) 

Classical CLT 


F (X) = % (X) = & 3qi -i (X) stable distribution 

1+91 

[with same (72 9 -i of /(X)] 
f 9{X) if |X|<X C ( 9 ,2); 

I C q ,2l\X\ 2 '^ if|X|»X c (^,2) 
for <7 > 1, with lim^i X c (q,2) = °° 






F (X) = Jzf^a (X) stable distribution 
[with same \X\ — ► °° behaviour of /(X)] 


Olq- 1 -» 00 

(0 < a < 2) 


F (X) = L a (X) stable distribution 
[with same \X — > °° behaviour of f(X)] 

( Sf(X) if |X|<X c (l,a); 
L a (X) ~ ^ 

I C 1:a /|X| a+1 if |X|»X C (?,2) 
with lim a ^2^t (1 5 Of) =°° 

Levy — Gnedenko CLT 


( 2 (l- ? )-o(3-«) 

sf2(.-,)-„ ( , + ,) a (x)~c„/|x| 

2(l- ? )-a(3-«5' a 

intermediate regime 

[x c (1) (4,a)«|X|<X e (2) ($,a)l; 

a+1 

^ a ,- r 3 2 (X)~C^/|X|'+«( 9 -D 
distant regime 

[|x|»xi 2) ( 9 ,o)l 

V L J 



considering q + -correlated variables (Eq. (12)) sharing the same distribution, and pre- 
senting a finite (2q — l)-variance, it is possible to observe the existence of a line 
of attraction, ^ q (X), (a = 2) in q — a space. We have also dealt with a ^-generalisation 
of the a-stable Levy distribution, Sf q , a (X) = [%,a (*)] (X). Since 3f q ,a(X) is 
stable, the convolution of such distributions, assuming they have an infinite (2q— 1)- 
variance, converges towards a J^f 9;Ct (X) distribution, which for large is equiva- 
lent to a &q L: 2 (X) distribution with qi = 2g ^~" +3 . Removing the restriction of infinite 
(2q— l)-variance, the convolution of £?q,a (X) distributions, asymptotically equivalent 

to <g q% a (X) ~ \X\~ (|X| -> oo and q x = 1+2a ^f a {q ~ l) ), leads to a % a (X) distribu- 
tion for which q = 2(i-g)-a(3_g) • These two asymptotic laws for ^f q ,a PO correspond 
to the existence of a crossover that we have depicted in Fig. 4. 

In both cases referred above, i.e., addition of random variables with finite and in- 
commensurable standard deviation, it is possible to define a triplet of parameters which 
characterises the attractor, the correlation, and the scaling rate, {P att , P C or,Pscl} U] rem- 
iniscent of the ^-triplet {q S em<lreh<lstat} conjectured in Ref. [15]. In that article, it was 
proposed that the values of {q S en,<lrel,<lstat} for say systems like long-range Hamilto- 
nian systems characterised by the interaction-decay exponent a and the dimension d 
would respect inequalities such as q re i,q s ta > 1 an d q se n < 1- Considering the convolu- 
tion of ^-independent random variables associated with £f 9l (X), it has been shown that 

the triplet of parameters {Patt,Pcor,Pscl} corresponds in fact to or simply 

{ik-iilkilk+i} following mapping relations. In this way, we can replace P, that stands 
for parameter, by q used to represent entropic indices. Hence, the triplet of parameters 
in actualy leading to the ^-triplet {q att , q C om<lscal}- Establishing a bridge between the 
triplet porposed in Ref. [15] and the triplet we have discussed, we argue that q^i, the 
entropic index for the attractor, should equal q stat . Taking into account Eq. (10), we can 
write 

*t-i=2 . (28) 

Qk+1 

It is noteworhy that the ^-triplet conjecture [15] was first observed by NASA using data 
from Voyager 1 related to the solar wind at the distant heliosphere [16], and also in a 
paradigmatic complex system such as a financial market [17]. For the solar wind obser- 
vations, it has been inferred the relations q stat + \/q re i = 2, and q rel + l/q sen = 2, which 
are consistent, within experimental error, with the results obtained from the NASA data 
set. Again, for reasons presented above, q stat = qk-i- Using Eq. (28) we have q re i =qk+i, 
and q sen = qk+3- These two cases of correspondence between {q s tat,<lreh<lsen} and 
{<lk*,<lk+T<lk>} represent strong candidates for the description of the ^-triplets for con- 
servative and dissipative systems. 

In addition, let us also mention two transformations that appear quite often in prob- 
lems discussed within nonextensive concepts, namely q a {q) =2 — q and q m (q) = l/q. 
In fact, these transformations, additive and multiplicative dualities respectively, have 
shown to be at the basis of the relations between entropic indices. Explicitly, if we ap- 
ply both transformations n times in a row, we obtain [q a qm] n (?) = f^n^fj- Looking 
to Eq. (26), we notice that sequences with a = 2 and n = 0, ±2, ±4, . . ., or a = 1 and 



n = 0, ±1, ±2, . . . coincide with [q a q m ] n - This fact is quite remarkable since it reveals a 
connection between the sequences that emerge from the ^-generalised central limit theo- 
rems and the dualities presented here above. The physical interpretation of these as well 
as other relations between entropic indices constitutes an interesting open challenge. 

Finally, let us mention the fact that the Central Limit Theorem results appear to also be 
applicable to the sum of deterministic variables [18]. It has been proved that, when the 
maximum Lyapunov exponent is positive, the sum of deterministic variables obtained 
from some dynamical process leads to a Gaussian distribution. Recently, studies on 
the sum of deterministic variables obtained from dissipative and conservative systems 
have been made. For the logistic map [19] (dissipative system) at the edge of the chaos 
(vanishing Lyapunov exponent) it has been verified that the intermediate part and the tail 
of the distribution are consistent with a ^-Gaussian distribution. Concerning conservative 
systems, studies on the Hamiltonian Mean Field model [20] at its metastable state have 
shown the emergence of non-Gaussian attracting probability density functions when the 
sum of velocities is performed. These resulting distributions have been numerically quite 
well approached by g-Gaussians over its whole range of values. Further analysis of both 
conservative and dissipative systems might clarify the emergence of a new generalisation 
of the central limit theorem, but for deterministic variables. 
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